Hypergeometric Functions and Carlitz Differential 
Equations over Function Fields 



Anatoly N. Kochubei 

Institute of Mathematics, 
National Academy of Sciences of Ukraine, 
Tereshchenkivska 3, Kiev, 01601 Ukraine 
E-mail: kochubei@i. com.ua 



1 



Abstract 



The paper is a survey of recent results in analysis of additive functions over function 
fields motivated by applications to various classes of special functions including Thakur's 
hypergeometric function. We consider basic notions and results of calculus, analytic 
theory of differential equations with Carlitz derivatives (including a counterpart of regular 
singularity), umbral calculus, holonomic modules over the Weyl-Carlitz ring. 
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1 INTRODUCTION 



Let K be the field of formal Laurent series t = with coefficients £j from the Galois 

j=N 

field F q , ^tv 7^ if t ^ 0, q = p v , v G Z + , where p is a prime number. It is well known that 
any non-discrete locally compact field of characteristic p is isomorphic to such K. The absolute 
value on K is given by \t\ = q~ N , |0| = 0. The ring of integers O = {t G K : \t\ < 1} is 
compact in the topology corresponding to the metric dist(£, s) = \t — s\. The absolute value | • | 
can be extended in a unique way onto the completion K c of an algebraic closure of K. 

Analysis over K and K c initiated by Carlitz 5 j and developed subsequently by Wagner, 
Goss, Thakur, the author, and others (see the bibliography in OEH]) is very different from the 
classical calculus. The new features begin with an appropriate version of the factorial invented 
by Carlitz - since the usual factorial i\, seen as an element of K, vanishes for i > p, Carlitz 
introduced the new one as 

A = -l] q ... [1] 9< ~\ [i] = x«'-x(i> 1), D = 1. (1.1) 

An important feature is the availability of many non-trivial F g -linear functions, that is such 
functions / defined on F 9 -subspaces K C K that 

f(ti + h) = fih) + f(t 2 ), f(at) = af(t), 

for any t,ti,t 2 G K , a G ¥ g . Such are, for example, polynomials and power series of the form 
^Ofct' 3 '', in particular, the Carlitz exponential 

oo q n 

n=o n 

and its composition inverse, the Carlitz logarithm 

OO n 

n=0 n 

where L n = [n][n — 1] • • • [1] (n > 1), Lq = 1. The notion of the Carlitz exponential obtained 
a wide generalization in the theory of Drinfeld modules (see [39 ). On the other hand, 
in various problems going beyond the class of F 9 -linear functions, an extended version of the 
Carlitz factorial (and its Gamma function interpolations) is used, so that D n can be seen as 
"an F g -linear part" of the full factorial; see [TH 13*9"] and references therein for the details. 

Among other special classes of F 9 -linear functions there are various polynomial systems 
(see below), an analog of the Bessel functions and Thakur's hypergeometric function 

[3*7| 13*51 13*9*1 ■ The latter is defined as follows. 

For n G Z + , a G Z, denote 

if a > 1; 

iia<0,n< -a; (1.4) 
if a < 0, n > —a. 
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Then, for dj, b{ G Z, such that the series below makes sense, we set 

(ai)„ • • • (a r ), 



,F s (a l5 . . . , a r ; fc, . . . , b s ; z) = V ^ n ^ (1-5) 

^ (»i)n • • • {b s ) n D n 



n=0 



Thakur [371 EE1 EH] has carried out a thorough investigation of the functions (1.5) and 
obtained analogs of many properties known for the classical situation. In particular, he found 
an analog of the hypergeometric differential equation. Its main ingredients are the difference 
operator 

(Au)(t) = u(xt) - xu{t) 

(an inner derivation of composition rings of F^-linear polynomials or more general F^-linear 
functions) introduced by Carlitz [5], the nonlinear (F^-linear) operator d — V o A, and the 
Fg-linear Frobenius operator tu = u q . For example, the function y = 2-^1(0, b; c; z) is a solution 
of the equation 

(A - [-o])(A - [-b\)y = d(A - [1 - c\)y. (1.6) 

Here we touch only a part of Thakur's results (he considered also hypergeometric functions 
corresponding to other places of ¥ q (x), a version of (1.5) with parameters from K and its 
extensions etc). 

The Carlitz exponential ec satisfies a much simpler equation of the same kind: 

de c = e a , (1.7) 

so that the operator d may be seen as an analog of the derivative. The operator r is an analog 

of the multiplication by t, so that A is the counterpart of t—. 

dt 

The same operators appear in the positive characteristic analogs of the canonical commu- 
tation relations of quantum mechanics fUEI]. in the analog of the Schrodinger representation 
we consider, on the Banach space Cq(0, K c ) of continuous F^-linear functions on O, with values 
from K c (with the supremum norm), the "creation and annihilation operators" 

a + = t — I, a~ = d 

(I is the identity operator). Then 

a-a + - a + a- = [l] 1/q I, (1.8) 

the operator a + a~ possesses the orthonormal (in the non- Archimedean sense |34j ) eigenbasis 
{fih 

(a + a~)/, = z = 0,l,2,...; (1.9) 
a + and a~ act upon the basis as follows: 

a + /i_! = [i]f u a~fi = i > 1; a~f = 0. (1.10) 

Here {f{\ is the sequence of normalized Carlitz polynomials 

A ri n ( s - m ) t^ 1 )' /o(s)= s , (1.11) 

mGF g [x] 
degm<i 



which forms an orthonormal basis in Cq(0, K c ). The spectrum of the "number operator" a + a~ 
is the set of elements [i], so that even this notation (proposed by Carlitz in 1935) becomes 
parallel to the usual quantum mechanical situation. 

An analog of the Bargmann-Fock representation is obtained if we consider the operators of 
almost the same form, 

a + = t, aT = d, 

but on the Banach space H of power series u(t) = Y) a n ^- with a n G K c , a n — > as n — > oo. 

— n n 

n=0 

These new operators satisfy the same relations (1.8)- (1.10), but this time instead of the Carlitz 
polynomials f n we get the eigenf unctions f n = 

The above results motivated the author to begin to develop analysis and theory of differential 
equations for F^-linear functions over K and K c , that is for the case which can be seen as a 
concentrated expression of features specific for the analysis in positive characteristic. This 
paper is a brief survey of some achievements in this direction. In particular, we consider the 
counterparts of the basic notions of calculus, analytic theory of differential equations (in the 
regular case and the case of regular singularity), their applications to some special functions, 
like the power function, logarithm and polylogarithms, Thakur's hypergeometric function etc. 
An umbral calculus and a theory of holonomic modules are initiated for this case. Like in the 
classical situation (see [7j), it is shown that some basic objects of the function field arithmetic 
generate holonomic modules. 

Note that some of the results can be easily extended to the case where the base field is a 
completion of F g (x) with respect to a finite place determined by an irreducible polynomial n G 
¥ q [x] (the field K corresponds to ir(x) = x); for some details see [23- The situation is different 
for the "infinite" place widely used in function field arithmetic (see [39J). In this case some of 
the basic objects behave in a quite different way - absolute values of the Carlitz factorials D n 
grow, as n —>■ oo, the Carlitz exponential is an entire function, the Carlitz polynomials do not 
form an orthonormal basis etc. A thorough investigation of properties of the Carlitz differential 
equations for this situation has not been carried out so far. 



2 Calculus 

2.1. Higher Carlitz operators A^ are introduced recursively, 

(AW«) (t) = A^u(xt) - x^A^uit), n>2. (2.1) 

For n = 1, the formula (2.1) coincides with the definition of A = A^\ if we set A^ ^ = I. 

The first application of these operators is the reconstruction formula ^H] for the coefficients 
a n of a power series u G H. Note that the classical formula does not make sense here because 
it contains the expression u^ n '(t)/n\ where both the numerator and denominator vanish. 



Theorem 2.1. If u G H , then 



A^u(t) 
lim jA±, n = 0,l,2, 
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For a continuous non-holomorphic F q -linear function u the behaviour of the functions 

® k u(t) =t- qk A^u{t), tGO\{0}, 

near the origin measures the smoothness of u. We say that u E Cq +1 (0,K c ) if D k u can 
be extended to a continuous function on O. This includes the case (k = 0) of different iable 
functions. 

The next theorem proved in ^H] gives a characterization of the above smoothness in terms 
of coefficients of the Fourier-Carlitz expansion. It includes, as a particular case (k = 0), the 
characterization of differentiable F g -linear functions obtained by Wagner [32|. 

oo 

Theorem 2.2. A function u — YJ c nfn £ Co(0, K c ) belongs to Cq +1 (0, K c ) if and only if 

n=0 

k 

Q nq \c n \ — * for n — > oo. 

In this case 

sup|D fc M(t)| = supg (n - fc)<?fe |c n |. 

teO n>k 

For a generalization to some classes of not necessarily F g -linear functions see |45j . 
Similarly ^Hl, a function u is analytic on the ball O (that is, u(t) = Yl a it q \ a i ~^ 0) if 

and only if q"- 1 !^! — > 0, as n — > oo. A more refined result by Yang jUj, useful in many 
applications, which makes it possible to find an exact domain of analyticity, is as follows (again 
we consider only Fg-linear functions while in jH] a more general class is studied). 

oo 

Theorem 2.3 (Yang). A function u = ^2 c n f n G C (O,K c ) is locally analytic if and only if 

n=0 

7 = liminf {-g- n logJc n |} > 0, (2.2) 

n^oo 

and if (2.2) holds, then u is analytic on any ball of the radius q~ l , 

I = max(0, [— (log(g - 1) + log 7)/ logg] + 1). 

2.2. Viewing d as a kind of a derivative, it is natural to introduce an antiderivative S 
setting Sf = u where u is a solution of the equation du = f, with the normalization u(l) = 0. 
It is easy to find Sf explicitly if / is given by its Fourier-Carlitz expansion (see [T5]). 

Next, we introduce a Volkenborn-type integral of a function / 6 Cq(0,K c ) (see for a 
similar integration theory over Z p ) setting 

[f(t)dt't i lim^l = (Sfy(0). 
o 

The integral is a F g -linear continuous functional on Cq(0, K c ), 

J cf(t)dt = c q J f(t)dt, ceK c , 
o o 
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possessing the following "invariance" property (related, in contrast to the case of Z p , to the 
multiplicative structure): 

r f(xt)dt = x [ f(t)dt-r(i). 



o o 

Calculating the integrals of some important functions we obtain new relations between them. 
In addition to the Carlitz exponential ec and the Carlitz polynomials f n (see (1.2) and (1.11)), 
we mention the Carlitz module function 

oo 

C s (z) =^2f i (s)z" i , seO,\z\<l. (2.3) 

i=0 

Note that if s G F g [s], then only the terms with i < degs are different from zero in (2.3). 
We have 

t<i n dt = -- r =*, n = 0,1,2, 

[n + 1\ 



[n -t- ij 

o 

\"+ 1 

n — f, 

5 5 5***5 



f "I " + i 

/ fn(t)dt= K —^ , 71 = 0,1,2,... 

O 

J C s {z) ds = \og c (z) — z, z £ K, \z\ < 1; 
o 

y e c (st) ds = t- ec{t), t <E K, \t\ < I. 



o 



For the proofs see 



3 Differential Equations for F g -Linear Functions 

3.1. Let us consider function field analogs of linear differential equations with holomorphic or 
polynomial coefficients. Note that in our situation the meaning of a polynomial coefficient is 
not a usual multiplication by a polynomial, but the action of a polynomial in the operator r. 
We begin with the regular case and consider an equation (actually, a system) 

dy{t) = P(r)y{t) + f(t) (3.1) 

where for each z 6 (i^ c ) m , t G K, 

P(r)z = J2^ k , /(*) = (3-2) 

k=0 j=0 3 

TXk are m x m matrices with elements from K c , (pj G (f^ c ) m , and it is assumed that the series 
(3.2) have positive radii of convergence. The action of the operator r upon a vector or a matrix 

is defined component- wise, so that z qk = \z\ , . . . , J for z — (zi, . . . , z m ). 
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We seek a F g -linear solution of (3.1) on some neighbourhood of the origin, of the form 



oo 



v(*) = X>|U y^{K c ) m , (3.3) 

1=0 

where yo is a given element, so that the "initial" condition for our situation is 

\imr l y(t) = y . (3.4) 

The next theorem, proved in jTHJ , is the function field analog of the Cauchy theorem from 
the classical analytic theory of differential equations. 

Theorem 3.1. For any y G (K c ^ m the equation (3.1) has a unique local solution of the form 
(3.3), which satisfies (3.4), with the series having a positive radius of convergence. 

Thus, regular equations with Carlitz derivatives behave more or less as their classical coun- 
terparts. The situation is different for singular equations. Let us consider scalar equations of 
arbitrary order 

m 

J2Mr)d'u = f (3.5) 

j=0 

oo £Q" 

where f{t)= fn — > A?'( r ) are power series having (as well as the one for /) positive radii 

n=o D n 
of convergence. 

We investigate formal solutions of (3.5), of the form 

u(t) = 22u n —, u n eK c . (3.6) 

n=0 n 

oo 

One can apply an operator series A(r) = a k rk (even without assuming its convergence) to 

k=0 

a formal series (3.6), setting 



r k u(t) = [n + if" ...[n + k]^—, k > 1, 

and 



Dn+k 



oo q l 

A{r)u{t) = E rf E ^ u il n + •■> + *] 

1=0 1 n+k=l 

where the factor [n + l] qk 1 . . . [n + k] is omitted for k = 0. These formal manipulations are 
based on the identity 

U 



Using also the relation 



d 



t q \ t q 



DiJ A-i' 

now we can give a meaning to the notion of a formal solution of the equation (3.5). 
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Theorem 3.2. Let u(t) be a formal solution (3.6) of the equation (3.5), where the series for 
Aj(r)z, z G K c , and f(t), have positive radii of convergence. Then the series (3.6) has a 
positive radius of convergence. 

This result (proved in is in a strong contrast to the classical theory. Note that in the 
p-adic case a similar phenomenon takes place for equations satisfying certain strong conditions 
upon zeros of indicial polynomials PjJ |H1 123 E2] ■ In our case such a behavior is proved for any 
equation, which resembles the (much simpler) case j2Zj of differential equations over a field of 
characteristics zero, whose residue field also has characteristic zero. 

3.2. The equations (3.1) and (3.5) behave like linear equations, though they are actually 
only Fq-linear. Theorem 3.1 can be extended [22| t° the case of strongly nonlinear equations 
(containing self-compositions y o y o ■ ■ ■ o y). 

On the other hand, it is natural to consider some equations of this kind in wider classes of 
Fg-linear functions resembling meromorphic functions of a complex variables. The set TZk of 
locally convergent F g -linear holomorphic functions forms a non- commutative ring with respect 
to the composition operation (the pointwise multiplication violates the F g -linearity) . The non- 
commutativity of TZk makes the algebraic structures related to Carlitz differential equations 
much more complicated compared to their classical counterparts. So far their understanding 
is only at its initial stage. It is known, however, that TZk can be imbedded into a skew field 
of F g -linear "meromorphic" series containing terms like t q (see j22] )■ A deep investigation of 
bi-infinite series of this kind convergent on the whole of K c has been carried out by Poonen 

A specific class of equations with solutions meromorphic in the above sense is the class of 
scalar Riccati-type equations 

dy(t) = \(y o y){t) + (P(r)y)(t) + R(t) (3.7) 

where A G K c , 

oo oo 

(P(r)y)(t) = £>/(*), R(t) = J2 r kt q \ 

k=l k=0 

Pk,Tk G K c (note that the right-hand side of (3.7) does not contain the linear term). The 
following theorem is proved in [22] • 

Theorem 3.3. IfO < |A| < q' 1/q2 , \p k \ < q~ 1/q2 , \r k \ < q- l/q2 for all k, then the equation (3.7) 
possesses solutions of the form 

oo 

y (t) = ct l ' q + J2 a nt gn , c,a n e K c , C + 0, 

n=0 

where the series converges on the open unit disk \t\ < 1. 
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4 Regular Singularity 



4.1. In analysis over C, a typical class of systems with regular singularity at the origin ( = 
over C consists of systems of the form 

(y'(C)= \B + ^A k Ay(() (4.1) 

where B, Aj are constant matrices, and the series converges on a neighbourhood of the origin. 
Such a system possesses a fundamental matrix solution of the form W(QC C where W(() is 
holomorphic on a neighbourhood of zero, C is a constant matrix, ( c = exp(ClogC) is defined 
by the obvious power series. Under some additional assumptions regarding the eigenvalues of 
the matrix B, one can take C = B. For similar results over C p see |llj . 

In order to investigate such a class of equations in the framework of Fq-linear analysis over K, 
one has to go beyond the class of locally analytic functions. Instead of power series expansions 
we can use the expansions in Carlitz polynomials on the compact ring O C K. The property 
of local analyticity, if it takes place, can be recovered with the use of Theorem 2.3. Note that 
our approach would fail if we consider equations over K c instead of K (our solutions may take 
their values from K c , but they are defined over subsets of K). In this sense our techniques are 
different from the ones developed for both the characteristic zero cases. 

We begin with the simplest model scalar equation 

rdu = Xu, A G K c , (4.2) 

whose solution may be seen as a function field counterpart of the power function t ^ t x . 

We look for a continuous F g -linear solution u(t, X) of the equation (4.2), with the "initial 
condition" u(l, X) = 1, in the form 

oo 

u(t) = X)c i / f (t) J teO, (4.3) 

i=0 

where c = 1. 

It is easy to see that the equation (4.2) has no continuous solutions if |A| > 1. If |A| < 1, 
then the solution u(t, A) is unique, continuous on O, and the coefficients from (4.3) have the 
form 

n-l 

Cn = H(X-\j})- 

j=0 

The function u(t, A) is analytic on O if and only if A = [j] for some j > 0; in this case 
u(t, A) = u(t, [j]) = t qJ . If A ^ [j] for any integer j > 0, then u(t, A) is locally analytic on O if 
and only if A = —x, and in that case u(t, —x) = for \t\ < q~ x . The relation 

u(t qm ,X) = u(t } X qm + [m\), teO, 

holds for all A, |A| < 1, and for all m — 0, 1, 2, For the proofs see [2*0] . 
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Similarly, if in (4.1) A = (A^) is is a m x m matrix with elements from K c , and we look 
for a matrix- valued solution of (4.1), then such a solution is given by the series (4.3) with the 
matrix coefficients 

Ci= jn(A-[j]/jje , «>1 

(I m is the unit matrix), if |A| = f max |A^| < 1. 

4.2. The analog, for our situation, of the system (4.1) is the system 

rdu - P(t)u = (4.4) 

oo k 

where P(t) is a matrix-valued analytic function, so that P(t)z = ^2 7TkZ q . We assume that 

k=0 

Kfc| < 7) 7 > 0, for all k, \n \ < 1. Denote by g(t) a solution of the equation rdg = n Q g. Let 
Ai, . . . , A m G K c be the eigenvalues of the matrix tt . 

Theorem 4.1. If 

Xi-\f^[k], i,j = l,...,m; k = l,2,..., (4.5) 
then the system (4-4) has a matrix solution 

oo 

u(t) = W(g{t)), W(s) = ^w k s q \ w = I m , 

where the series for W has a positive radius of convergence. 

The paper [20J contains, apart from the proof of Theorem 4.1, a discussion of some situations 
(the Euler type equations) where its conditions are violated, as well as of the meaning of the 
conditions (4.5). Here we only mention that in the scalar case m = 1 the condition (4.5) is 
equivalent to the assumption 7r ^ —x, so that it excludes the case where solutions of the 
equation rdg = 7i g has pathological properties. 

4.3. For the above equation, continuous solutions were found as Fourier-Carlitz expansions 

oo 

u(t)=Y,Cnfn(t), (4.6) 
n=0 

and we had to impose certain conditions upon coefficients of the equation, in order to guarantee 
the uniform convergence of the series on O (which is equivalent to the fact that c n — ► 0). 
However formally we could write the series (4.6) for the solutions without those conditions. 
Thus, it is natural to ask whether the corresponding series (4.6) converge at some points t G O. 
Note that (4.6) always makes sense for t G W q [x] (for each such t only a finite number of terms is 
different from zero). The question is whether the series converges on a wider set; if the answer 
is negative, such a formal solution is called strongly singular. 

The available results regarding strong singularity of solutions of some equations are based 
on the following general fact [2"Uj . 
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Theorem 4.2. If |cj| > p > /or a// i > io (where io is some natural number), then the 
function (4-6) is strongly singular. 

It follows from Theorem 4.2 that non-trivial formal solutions of the equation (4.2) with 
|A| > 1 are strongly singular. A more complicated example is provided by the equation 

(A - [-a]) (A - [-b])u = dAu, a,beZ, (4.7) 

for Thakur's hypergeometric function 2-^1(0, b; l;i). 

A holomorphic solution of (4.7) is given by an appropriate specialization of (1.5). Classically 
(over C), there exists the second solution with a logarithmic singularity. Here the situation is 
different. Looking for a solution of the form (4.6) we obtain a recursive relation 

(c& - c l+2 ) + c%\i + l]^ - c m ([.] + [i + 1] - [-a] - l-b}) 

-Ci([i]-[-a])([z]-[-6]) = 0, z = 0,l,2,.... (4.8) 

Taking arbitrary initial coefficients Co, C\ £ K c we obtain a solution u defined on ¥ g [x}. On 
each step we have to solve the equation 

z yi - z = v. (4.9) 

If |cf| < 1 and |cj + i| < 1, then in the equation for c i+2 we have \v\ < 1. 

It can be shown [20] that the equation (4.9) has a unique solution z £ K c , for which 
kol — 1^1) an d — 1 other solutions z, \z\ = 1. It is natural to call a solution generic if, starting 
from a certain step of finding the coefficients c n , we always take the most frequent option 
corresponding to a solution of (4.9) with \z\ = 1. Now Theorem 4.2 implies the following fact. 

Theorem 4.3. A generic solution of the equation (4-7) is strongly singular. 

Of course, in some special cases the recursion (4.8) can lead to more regular solutions, in 
particular, to the holomorphic solutions found by Thakur. 



5 Poly logarithms and a Zeta Function 

5.1. The Carlitz differential equations can be used for defining new special functions with 
interesting properties. Some examples are given in this section. 

An analog of the function — log(l — t) is defined via the equation 

(1 - r)du(t) = t, teKi,, (5.1) 

a counterpart of the classical equation (1 — t)u'(t) = 1. The next results are taken from [2*3] 
where the equation (5.1) is considered for an arbitrary finite place of ¥ q (x). 

Let li(t) be a F g -linear holomorphic solution of (5.1) with the zero initial condition (in the 
sense of (3.4)). Then it is easy to show that 

OO q n 
71=1 L J 
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and the series in (5.2) converges for \t\ < q^ 1 . 

Note that l\(t) is different from the well-known Carlitz logarithm log c (see (1.3)). Analogies 
motivating the introduction of special functions are not so unambiguous, and, for instance, from 
the composition ring viewpoint, log c is an analog of e _< , though in other respects it is a valuable 
analog of the logarithm. By the way, another possible analog of the logarithm is a continuous 
function u(t), \t\ w < 1, satisfying the equation Au(t) = t (an analog of tu'(t) = 1) and the 
condition u(l) = 0. In fact, u — T>\, the first hyperdifferential operator (the definition of T>\ is 
given in Sect. 5.2 below); see [To] . 

Now we consider continuous non-holomorphic extensions of l\. 

Theorem 5.1. The equation (5.1) has exactly q continuous solutions on O coinciding with 

oo 

(5.2) as \t\ < q~ l . These solutions have the expansions in the Carlitz polynomials u = Y2 c ifi 

i=0 

where C\ is an arbitrary solution of the equation c\ — C\ + 1 = 0, higher coefficients are found 
from the relation 

oo 

C n = 22 ( C n-ll n ~ i]) 9 ' , n > 2, 
3=0 

and the coefficient cq is determined by the relation 

oo 

co = ^(-^^f > 

i=l 

Below we denote by l\ an arbitrary fixed "branch" of extensions of (5.2). 
The polylogarithms l n (t) are defined recursively by the equations 

Al n = /„_!, n > 2, (5.3) 

which agree with the classical ones tl' n (t) = l n -i(t). Analytic F g -linear solutions of (5.3), such 
that t~ l l n (i) — > as t — > 0, are found easily by induction: 

OO gj 

W*) = Er^' (5 ' 4) 



Theorem 5.2. For each n > 2, there exists a unique continuous ¥ q -linear solution of the 
equation (5.3) coinciding for \t\ < q^ 1 with the polylogarithm (5.4). The solution is given by 

°° (n) 

the Carlitz expansion l n — ^2 c \ fi with 

i=0 

<C n q-*-\ C n >0, i>l, 

oo (n) 
i=l 




13 



5.2. Now that the above polylogarithms have been extended onto the disk {\t\ < 1}, we can 
interpret their values at t = 1 as "special values" of a kind of a zeta function. In order to define 
the latter, we introduce the operator A^ a \ a G O, a function field analog of the Hadamard 
fractional derivative (t-^) a from real analysis (see [33]). 

Denote by V k {t), k > 0, t G O, the sequence of hyperdifferentiations defined initially on 
monomials by the relations V (x n ) = x n , X>fc(l) = 0, k > 1, 



W 4 n-fc 



V k (x n ) = \ k Jx 

where it is assumed that (?) = for k > n. V k is extended onto ¥ q [x] by Fq-linearity, and 
then onto O by continuity [12] • The sequence {T>k} is an orthonormal basis of the space of 
continuous F^-linear functions on O jTHUHl- 

oo oo 

Let a G O, a = a n x n , a n G ¥ q . Denote a = (~ l) n a n x n . For an arbitrary continuous 

n=0 n=0 

Fq-linear function u on O we define its "fractional derivative" A^u at a point t G O by the 
formula 

oo 

(A (Q) n) (t) = ^(-l) fc D fc (a)n(x fc t). 

fc=0 

The function a i— > (A^'m) (t) is continuous and F^-linear. As a function of t, A^u is continuous 
if, for example, u is Holder continuous. 

Our understanding of A^ as a kind of a fractional derivative is justified by the following 
properties: 

A (xn) = A n , n = l,2,...; 
A (a) (A^u) (t) = (A (Q/3) n) (t), 

for any a, (3 G O. 

5.3. We define £(*), t G setting C(0) = 0, 

C(x- n ) = l n (l), n = l,2,..., 

and 

C(f) = (A^+^+-\) (1), n = l,2,... l 

if i = x~™(#o + + '••)) ^g- The function £ is a continuous F g -linear function on K x 
with values in K c . 

In particular, we have 

((x m ) = (A m+1 /i) (1), m = 0,1,2,.... 
The above definition is of course inspired by the classical polylogarithm relation 



*:r)E- = E 

az J L — ' n s z — ' n 

' n=l n=l 



s-1 



In contrast to Goss's zeta function defined on natural numbers and interpolated onto Z p (see 
[HUSH!), the above ( is purely an object of the characteristic p arithmetic. 
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Let us write some relations for special values of our £; for the details see 
As we saw, 

oo „j oo 

j = l 4 = 

Next, let us consider the double sequence A n>r 6 K, A n ^ = (— l) n-1 L n _i, 

A n , r = (-l) n + r L^ 1 F1F1~ r r r ^ 2 ' 



These elements appear as the coefficients of the expansion [32] of a hyperdifferentiation T> r in 
the normalized Carlitz polynomials, as well as in the expression of the operators A*™) from 
(2.1) via the iterations A r . Here we have the identity 

oo i 

c(o = E(- 1 ) i+lL r 1 EArC(a: r - B ) 

8=1 r=l 

which may be seen as a distant relative of Riemann's functional equation for the classical zeta. 

Finally, consider the coefficients q of the Carlitz expansion of l\ (see Theorem 5.1). They 
are expressed via zeta values: 



r=l 



By Theorem 5.1, for % > 2 we have 



3=0 

The series in (5.5) may be seen as an analog of ^ j~ z . This analogy becomes clearer if, for a 

j 

oo 

fixed z G K c , \z\ < 1, we consider the set S of all convergent power series ^ z q3 " corresponding 

n=l 

to sequences {j n } C N. Let us introduce the multiplication in S setting z q% z qJ = z q ' J and 
extending the operation distributively (for a similar construction in the framework of g-analysis 

in characteristic see [23 )■ Denoting by Y\ t ne product in S of elements indexed by prime 

p 

numbers we obtain in a standard way the identity 



n=0 



(the infinite product is understood as a limit of the partial products in the topology of K c ), an 
analog of the Euler product formula. It would be interesting to study the algebraic structure 
of S in detail. 
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6 Umbral Calculus 



6.1. Classical umbral calculus (HUEO] is a set of algebraic tools for obtaining, in a unified way, a 
rich variety of results regarding structure and properties of various polynomial sequences. There 
exists a lot of generalizations extending umbral methods to other classes of functions. However 
there is a restriction common to the whole literature on umbral calculus - the underlying field 
must be of zero characteristic. An attempt to mimic the characteristic zero procedures in the 
positive characteristic case [T2| revealed a number of pathological properties of the resulting 
structures. More importantly, these structures were not connected with the existing analysis 
in positive characteristic based on a completely different algebraic foundation. 

A version of umbral calculus inmplementing such a connection was developed by the author 
|2*T] , and we summarize it in this section. Its basic notion is motivated by the following identity 
for the non-normalized Carlitz polynomials = Diff. 

e l (st) = J2( t ) e n {t){e^ n {s)Y n (6.1) 

n=0 VV K 

where the "i^-binomial coefficients" I ] are defined as 



n, K 



n) K DnDtn 

Computing the absolute values of the Carlitz factorials directly from their definition (1.1), 
it is easy to show that 

1, < n < i. 



n J K 

In fact, ( ) e F»(x), and we can consider also other places of F„(x), that is other non- 

w* 

equivalent absolute values. It can be proved EH that I I belongs to the ring of integers for 

any finite place of ¥ q (x). 

We see the relation (6.1) as a function field counterpart of the classical binomial identity 
PH E?U] satisfied by many classical polynomials. Now, considering a sequence U{ of F g -linear 
polynomials with coefficients from K c , we call it a sequence of K -binomial type if degUj = q l 
and for all i = 0, 1, 2, . . . 

Ui( st )=J2( l ) u n (t) { Ui - n (s)} qn , s,teK. (6.2) 

n=0 K 

As in the conventional umbral calculus, the dual notion is that of a delta operator. However, 
in contrast to the classical situation, here the delta operators are only F^-linear, not linear. 

Denote by p\ the operator of multiplicative shift, (p\u)(t) = u(Xt). We call a linear operator 
T, on the i^ c -vector space K c {t} of all F g -linear polynomials, invariant if it commutes with p\ 
for each A G K. 



16 



A F q -linear operator 8 = r _1 <5o, where Sq is a linear invariant operator on K c {t}, is called 
a delta operator if 8o(t) = and 80(f) 7^ for deg / > 1. A sequence {P n }o° of Fg-linear 
polynomials is called a basic sequence corresponding to a delta operator 5 = t~ 1 8o, if degP n = 
q n , P (l) = 1, P n (l) = forn > 1, 

8P = 0, <5P„ = [nJ^Pn-x, 71 > 1, (6.3) 

or, equivalently, 

5 Po = 0, 5 P n = WCi, 7i > 1. (6.4) 

It is clear that d = r _1 A is a delta operator. It follows from well-known identities for the 
Carlitz polynomials [T3] (see also (1.10)) that the sequence {e^} is basic with respect to the 
operator d. 

Theorem 6.1. For any delta operator 8 = t _1 5 , there exists a unique basic sequence {P n }, 
which is a sequence of K -binomial type. Conversely, given a sequence {P n } of K -binomial type, 
define the action of 8 on P n by the relations (6.4), extend it onto K c {t} by linearity and set 
8 = r _1 5 . Then 5 is a delta operator, and {P n } is the corresponding basic sequence. 

The analogs of the higher Carlitz difference operators (2.1) in the present general context 
are the operators 8$ = r l 5 l . The identity 

^ = ^-Pf-l (6-5) 

holds for any I < j. If / is a F^-linear polynomial, deg/ < q n , then a generalized Taylor 
formula 

" (*?>/) 00 

f( st ) = I2 d m (6 - 6) 

holds for any s,t G K. For the Carlitz polynomials e^, the formulas (6.5) and (6.6) are well 
known [13J. It is important that, in contrast to the classical umbral calculus, the linear operators 
involved in (6.6) are not powers of a single linear operator. 

Any linear invariant operator T on K c {t} admits a representation 

rr jt(0 ( TP l)i l ) (a *7\ 

T = Zs (Tl °°> ° l = — n — • 



1=0 



The infinite series in (6.7) becomes actually a finite sum if both sides of (6.7) are applied to any 
Fg-linear polynomial. Conversely, any such series defines a linear invariant operator on K c {t}. 

Let us consider the case where 5 = d, so that 8^ = A^. The next result leads to new delta 
operators and basic sequences. 

Theorem 6.2. The operator 9 = t~ 1 6q, where 

00 

6>o = 5>AW 
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is a delta operator if and only if 



n 

def CTJ 



S n = 2^ — L_ ^ /or a// n = 1, 2, . . . . (6.1 
i=i D n _i 



Example 1. Let o~\ = 1 for all / > 1, that is 

oo 

e = J2 A{l) - ( 6 - 9 ) 



1=1 



Estimates of \D n \ which follow directly from (1.1) show that \S n \ = q i- 1 , so that (6.8) is 
satisfied. Comparing (6.9) with a classical formula from [3T] we may see the polynomials P n 
for this case as analogs of the Laguerre polynomials. 

Example 2. Let o\ = — . For this case it can be shown [21] that S n = D^, n = 

Li 

1,2,.. .; e {t qJ ) = t qj for all j > 1 (of course, o (t) = 0), and P (t) = t, P n {t) = D n (V" - t 9 "" 1 ) 
for n > 1. 

6.2. As in the p-adic case |IUJ |JTJ , ^ e um ^ ra ^ ca l cums can De used for constructing new 
orthonormal bases in Cq(0, K c ). 

Let {P n } be the basic sequence corresponding to a delta operator 5 = r _1 <5 , 

oo 

5o = J>A«. (6.10) 

i=i 

Pn 

The sequence Q n = — — , n = 0,1,2,..., called the normalized basic sequence, satisfies the 
identity 

i 

Qi(st) = ^2Q n (t) {Q t - n (s)} qn , 

n=0 

another form of the .ff -binomial property. Though it resembles its classical counterpart, the 
presence of the Frobenius powers is a feature specific for the case of a positive characteristic. 

Theorem 6.3. If\o~x\ = 1, \o~{\ < 1 fori > 2, then the sequence {Qn}^ ^ s an orthonormal basis 
of the space Cq(0, K c ) - for any f G Cq(0, K c ) there is a uniformly convergent expansion 

oo 

/(*) = J^»Qn(f), teo, 

n=0 

where ip n = (Sg 1 ^ fj (1), \ip n \ — > as n — > oo ; 

= SUp \lp n \. 
n>0 
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By Theorem 6.3, the Laguerre-type polynomial sequence from Example 1 is an orthonormal 
basis of Cq(0, K c ). The sequence from Example 2 does not satisfy the conditions of Theorem 
6.3. 

Note that the conditions of Theorem 6.3 imply that S n ^ for all n, so that the series 
(6.10) considered in Theorem 6.3 always correspond to delta operators. 

In |2T] recursive formulas and generating functions for normalized basic sequences are also 
given. 

7 The Weyl-Carlitz Ring and Holonomic Modules 

7.1. The theory of holonomic modules over the Weyl algebra and more general algebras of 
differential or g-difference operators is becoming increasingly important, both as a crucial part 
of the general theory of D-modules and in view of various applications (see, for example, jlj [Tj 
32j). Usually, the holonomic property of the module corresponding to a system of differential 
equations is a sign of its "regular" behavior. Most of the classical special functions are associated 
(see [7j) with holonomic modules, which helps to investigate their properties. 

It is clear from the above results that in the positive characteristic case a natural counterpart 
of the Weyl algebra is, for the case of a single variable, the ring 2li generated by r, d, and scalars 
from K c , with the relations 

dr - rd = [l} 1/g , rX = X q r, dX = X 1/q d (A G K c ). (7.1) 

The ring consists of finite sums 

<> X>,7'r/<. \ {j £ K c , (7.2) 

and the representation of an element in the form (7.2) is unique. 

Basic algebraic properties of 2ti ^3 HI] are similar to those of the Weyl algebra in charac- 
teristic and quite different from the case of the algebra of usual differential operators over a 
field of positive characteristic [28 . 

The ring 2li is left and right Noetherian, without zero divisors. 2li possesses no non-trivial 
two-sided ideals stable with respect to the mapping 

The centre of 21 1 is described explicitly in [3J; it contains countably many elements (this corrects 
an erroneous statement from jH]). In fact, 21 1 belongs to the class of generalized Weyl algebras 
j2]. A well-developed theory available for them enabled Bavula |3 j to classify ideals in 2t 1; as 
well as all simple modules over Six. 

A generalization of 2li to the case of several variables is not straightforward because the 
Carlitz derivatives d s and dt do not commute on a monomial f{s,t) = s qm t q, \ if m ^ n. 
Moreover, if m > n, then d™f is not a polynomial, nor even a holomorphic function in t (since 
the action of d is not linear and involves taking the q-th root). 
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A reasonable generalization is inspired by Zeilberger's idea (see [7j) to study holonomic 
properties of sequences of functions making a transform with respect to the discrete variables, 
which reduces the continuous-discrete case to the purely continuous one (simultaneously in 
all the variables). In our situation, if {Pfc(s)} is a sequence of F^-linear polynomials with 

degPfc < q k , we set 

oo 

/(*,*) = IC^OO***, (7-3) 

k=0 

and d s is well-defined. In the variable t, we consider not d t but the linear operator A t . The 
latter does not commute with d s either, but satisfies the commutation relations 

d s At ~ A t d s = [l]V*d s , A t r - rA t = [l]r, 

so that the resulting ring 2t 2 resembles a universal enveloping algebra of a solvable Lie algebra. 
More generally, denote by T n +\ the set of all germs of functions of the form 

oo oo min(fci,...,fc„) 

f(s, /, /„) Yl E ^ M ,.,k n s qm tf . . . C (7.4) 

fel=0 k„=0 m=0 

where a m ,fci,...,fc n G K c are such that all the series are convergent on some neighbourhoods of 
the origin. We do not exclude the case n = where T\ will mean the set of all F^-linear power 
series ^ a m s gm convergent on a neighbourhood of the origin. J- n+ \ will denote the set of all 

m 

polynomials from T n +\, that is the series (7.4) in which only a finite number of coefficients is 
different from zero. 

The ring 2l n+ i is generated by the operators r, d s , A tl , . . . A tn on JF n+1 , and the operators 
of multiplication by scalars from K c . To simplify the notation, we write Aj instead of 
and identify a scalar A G K c with the operator of multiplication by A. The operators Aj are 
i^ c -linear, so that 

AjX = XAj, A G K c , (7.5) 

while the operators r, d s satisfy the commutation relations (7.1). In the action of each operator 
d s ,Aj (acting in a single variable), other variables are treated as scalars. The operator r acts 
simultaneously on all the variables and coefficients. We have the relations involving A^-: 

A jT - rAj = [l]r, d s Aj - Ajd s = [l) 1/q d s , j = 1, . . . , n. (7.6) 

Using the commutation relations (7.1), (7.5), and (7.6), we can write any a G 2l n+ i, in a 
unique way, as a finite sum 

" ;,vVAV...A;;. (7.7) 

Let us introduce a filtration in (an analog of the Bernstein filtration) denoting by T u , 
v G Z + , the K c -vector space of operators (7.7) with max{/ + // + %\ + ■ ■ • + i n } < v where the 
maximum is taken over all the terms of (7.7). Then 2t n+ i is a left and right Noetherian filtered 
ring. 
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In a standard way (see [TU]) we define filtered left modules over 2t n+ i. All the basic notions 
regarding a filtered module M (like those of the graded module gr(M), dimension d(M), mul- 
tiplicity m(M), good filtration etc) are introduced just as their counterparts in the theory of 
modules over the Weyl algebra. 

If we consider 2L+ left module over itself, then 

d(On+i) = n + 2, m(2t„, +1 ) = 1. (7.8) 

For any finitely generated left 2l n+ i-module M, we have d(M) < n + 2. By (7.8), this bound 
cannot be improved in general. However, if J is a non-zero left ideal in 2l n+ i, then 

d(SUi//)<n + l. (7.9) 

For the module T n +\ of F^-linear polynomials (7.4), we have 

d (fn+i) =n + l, m (pn+i) = n\ 
The proofs of all these results, as well as the ones given in this section below, can be found in 

El- 

It is natural to call an -module M holonomic if d(M) = n + 1. Thus, !F n +i is an 
example of a holonomic module. 

The next theorem demonstrates, already for the case of 2li-modules, a sharp difference from 
the case of modules over the Weyl algebras. In particular, we see that an analog of the Bernstein 
inequality (see [TU]) does not hold here without some additional assumptions. 



Theorem 7.1. (i) For any k = 1,2,..., there exists such a nontrivial Qlx-module M that 
dimM = k (dim means the dimension over K c ), that is d(M) = 0. 

(ii) Let M be a finitely generated %\-module with a good filtration. Suppose that there exists a 
"vacuum vector" v G M, such that d s v = and r m (v) ^ for all m = 0, 1, 2, . . .. Then 
d(M) > 1. 

7.3. Let us consider the case of holonomic submodules of the 2l n+ i-module .F n +i, consisting 
of F g -linear functions (7.4) polynomial in s and holomorphic near the origin in t±, . . . , t n . 
Let ^ / G F n+ll 

I f = We SUj : ifU) = 0} . 

If is a left ideal in 2l n+ i. The left 2t n+ i-module Mf = ^i n+1 /If is isomorphic to the submodule 
2t„ + i/ C J-'n+i - an element <f(f) G 2l n+ i/ corresponds to the class of G 2t n +i in Mf. A 
natural good filtration in Mf is induced from that in 2l n+ i. 

As we know (see (7.9)), if If ^ {0}, then d(Mf) < n+1. We call a function / holonomic if 
the module Mf is holonomic, that is d(Mf) = n + 1. The condition If ^ {0} means that / is 
a solution of a non-trivial "differential equation" ip(f) = 0, (p G 2l n +i- The case n = is quite 
simple. 

Theorem 7.2. If a non-zero function f G T\ satisfies an equation <p(f) = ; ^ cp G 2li, i/ien 
/ zs holonomic. 
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In particular, any F g -linear polynomial of s is holonomic, since it is annihilated by d™, with 
a sufficiently large m. 

If n > 0, the situation is more complicated. We call the module Mf (and the corresponding 
function /) degenerate if D(Mf) < n + 1 (by the Bernstein inequality, there is no degeneracy 
phenomena for modules over the complex Weyl algebra). The simplest example of a degenerate 
function (for n = 1) is f(s,ti) = g(sti) G Ti where the function g belongs to T\ and satisfies 
an equation <p(g) — 0, <p G 2li. It can be shown that d{Mf) = 1. 

In order to exclude the degenerate case, we introduce the notion of a non-sparse function. 

A function / G J- n +i of the form (7.4) is called non-sparse if there exists such a sequence 
mi —> oo that, for any I, there exist sequences kf\ k%\ ■ ■ ■ , kn > mi (depending on I), such 
that ku —> oo as i —>■ oo (z/ = 1, . . . , n), and a (o 7^ 0. 

Theorem 7.3. J/ a function f is non-sparse, then d(Mf) > n + 1. //, m addition, f satisfies 
an equation (p(f) = 0, 7^ cp G 2l n+ i, i/ien / zs holonomic. 

7.4 We use Theorem 7.3 to prove that the functions (7.4) obtained via the sequence-to- 
function transform (7.3) or its multi-index generalizations, from some well-known sequences of 
polynomials over K are holonomic. In all the cases below the non-sparseness is evident, and 
we have only to prove that the corresponding function satisfies a non-trivial Carlitz differential 
equation. 

a) The Carlitz polynomials. The transform (7.3) of the sequence {fk} is the Carlitz module 
function C s (t); see (2.3). It is easy to check that d s C s (t) = C s (t). Therefore the Carlitz module 
function is holonomic, jointly in both its variables. 

b) Thakur's hypergeometric polynomials. We consider the polynomial case of Thakur's 
hypergeometric function (1.5), that is 

rn / u 7 \ \^ [~ a l)m ■ ■ ■ {— a l)m q m /» 1f >\ 

iF x (-a u . . . , -a { -61, . . . , -6 A z) = y —t- z 9 (7.10) 

^ (-Ol)m • • • {-0\) m D m 

where a\, . . . ,ai,b\, . . . ,b\ G Z + . It is seen from (1.4) that the terms in (7.10), which make sense 
and do not vanish, are those with m < min(ai, . . . ,ai,bi, . . . ,b\). Let the function / G Fi+x+i 
be given by 

f(s,t 1 ,...,ti,u 1 ,...,u x ) 

= Yl ■ ■ ■ Yl Yl ■ ■ ■ Yl i F ^(~ k ^ • • • ' -hi _z/ i' • • • ' _z/ ^; s) 

fc 1= ki=0vi=0 u x =0 



xtf ...tfuf ...u^ 



It is known (|^|, Sect. 6.5) that 
d s iF x (-ki, -ku -vi, -v X ] s) 



F x (-h + 1, . . . , -h + 1; -U! + 1, . . . , -u x + 1; s) (7.11) 
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if all the parameters ki, . . . , ki, ux, . . . , v\ are different from zero. If at least one of them is 
equal to zero, then the left-hand side of (7.11) equals zero. This property implies the identity 
d s f = f, the same as that for the Carlitz module function. Thus, / is holonomic. 

f k\ 

c). K -binomial coefficients. It can be shown [21] that the .ff -binomial coefficients 



[see Sect. 6) satisfy the Pascal-type identity 



m J K \m-lj K \ m 7 K 
where < m < k and it is assumed that I 1=1, 1=0. 



1/k V & / K 



Consider a function / e JF 2 associated with the i^-binomial coefficients, that is 



The identity (7.12) implies the equation 

d a f(s,t) = A t f(s,t) + [l] 1 '«f(s,t) 
for the function (7.13). Therefore / is holonomic. 



mj K 



k\ (k- l\ q /Ife-lV „, , 

+ ( ) D^ 1 (7.12) 



oo k / r \ 
k=0 m=0 ^ ' K 
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